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Abstract. Let Xd(p) De the p monomer-dimer entropy on the d-dimcnsional integer lattice Z d , where p £ 
[0, 1] is the dimer density. We give upper and lower bounds for \d{p) in terms of expressions involving 
Ad_i(g). The upper bound is based on a conjecture claiming that the p monomer-dimer entropy of an 
infinite subset of Z d is bounded above by \d(p)- We compute the first three terms in the formal asymptotic 
expansion of \d(p) m powers of ^. We prove that the lower asymptotic matching conjecture is satisfied for 
A d (p). 

Keywords and phrases: monomer-dimer entropies: asymptotic expansion, recursive inequalities and a 
conjecture. 

1. Introduction 

The first aim of this paper is to discuss an asymptotic expansion of the monomer-dimer p-entropy, denoted 
by X d (p), where p £ [0, 1] is the density of dimers, on the integer d-dimensional lattice Z d . The study of 
the existence of this entropy, some of its properties and its estimates, was initiated in a series of papers by 
Hammersley and his collaborators []j3 [201 HD 122] ■ It was shown by Mine [53] that the d-dimensional dimer 
entropy A^ (1) satisfies 

1, , 1 . , log (2d)! 1 . rt 1 log27r(2d) 1 . . 

- log (2d) - - < X d (1) < < - log (2d) ~ - + — i — - + jr. 1.1) 

2 6 v ' 2 ~ y J ~ id ~ 2 B v ' 2 id 48d 2 v ' 

The lower bound is implied by the proof of the van der Waerden conjecture [5] [6], or its weak form [15j . The 

first upper bound follows from the Bregman ineqality [2]. The last upper bound follows from a sharp form 

of the Stirling formula for (2d)! [T31 p 52]. In a series of papers [7] - [TT], the first author studied a possible 

asymptotic expansion of 

1 1 00 

A e2 (l)~-log(2d)--+£-|. (1.2) 

fc=i 

This is derived assuming that an argument employing a formal cluster expansion could be made rigorous. He 
was able to compute the values of c%, C2, C3, equal to |, A respectively. In this paper we show that the 
methods in [7] - [11] can be also applied to the asymptotic expansion of (p) for any p £ [0,1]. The lower 
and the upper bounds for Xd (p) given in [TT1 (5.3)] and [HI (7.5)] respectively yield 

~ (p log 2d - p logp - 2 (1 - p) log (1 - p) - p) < X d (p) < (1.3) 
\ (p^^-plogP-2(l-p)log(l-p)) < (1.4) 

i(plog2d-plogp-2(l-p) log (1 - p) - p) +p 2 ^ 2d) + -1^ 

Assuming valid the same formal procedure as in [7] or QT| , we show that 

1 °° ( ) 

Arf (p) ~ g (P log 2d - plogp - 2 (1 - p) log (1 - p) - p) + £ (1.5) 

fc=i 
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where each Ck (p) is a polynomial in p. This polynomial, Ck (p) , is built of powers p s with k < s < 2k ( as 
follows from the form of equation (|5.22[) below ) . 
We compute the polynomials C\ (p) , c-i (p) , C3 (p), 

ci (p) = y (1.6) 

c 2 (p) = (2 P 3 + 3p 4 ) /96 (1.7) 
c 3 (p) = (-5p 4 + 12p 5 + 8 P 6 ) /192. (1.8) 

While studying Ad (p) we made a heuristic observation 

A, (p) - A d _a (p (1 - i)) + (1 - p (1 - i)) log (1 - p (1 - I)) - £ log (£) 

+ (1 - log (1 - ^) - (1 - p) log (1 - p) , (1.9) 

for large d. We also substituted the ansatz 

\ d (p) = iplog (2d) - | logp - (1 - p) log (1 - p) - I + ^ + 5/ 9 + 2 Cp3 (1.10) 

into both sides of (|1.9[) . and found the value of the right side of the equation minus the left side of the 
equation is equal to 

4<- i+8 «>£-^<- 2+2f * + 4 +o (i) < lii > 

The correct values a — | and c = 2 are determined to achieve accuracy up to order ^3 (other choices would 
have terms with lower powers of p or -3). We view this result as both a measure of how good an approximation 
(|1.9I) is, and a verification that the series (| 1 . 5|) is correct (to a certain order). 
We show in this paper that (jl.9l) is actually a recursive inequality 

A, (p) > A<2_i (p (1 - I)) + (1 - p (1 - i)) log (1 -p (1 - i)) - £ log (£) 

+ (1 - &) (1 - ^) - (1 - P) log (1 - p) . (1.12) 
We also show the validity of some upper bounds on Ad (p) in terms of the function Xd—i (p)- 
A d (p) < 

(1 - I) log (1 - f) - f log § - (1 - u) log (!-«) + (!- u) A d _! (^^— ^ J (1.13) 



max 

«e[o,5] 



and 

A d (p) < 

max /(l-p+|)log(l-p+f) -|b g a-(l-p)log(l-p) + (l-§)A^i ff- 1 ^)) (1.14) 
«e[o,f] L - - ' V 1 _ 2 / J 

if the following conjecture holds. 

Conjecture 1.1. Let £ = {Si,i G N} fee an increasing sequence of finite sets in 7L d . We can define its p 
monomer- dimer entropy A5 (p) similar to the definition of Ad (p). We conjecture that A5 (p) < Ad (p). 

This conjecture will be precisely stated in Section [2] which contains some preliminary formalism. The lower 
bound is studied in Section [31 Section H] deals with the upper bounds, and goes a long way towards proving 
that the bound of (|1.14[) is always better than the bound of fll . 13[) . as we believe. The asymptotic expansion 
is developed in Section [5j Section [6] studies some implications of our inequalities, in particular proving the 
LAMC for the monomer-dimer problem on a rectangular lattice, see [16]. Section [7] studies numerical tests of 
our bounds and expansion for d = 1, 2, 3. Although our expansion was developed as an asymptotic expansion 
in 1 jd it has surprising validity for small d. 
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2. Preliminary Results 

Let Z d be the d-dimensional integer lattice in M. d . Two points x = (x\, . . . , x<j) , y = (j/i, • • • , Ud) £ Z d are 
called neighbors if ||x — y||i = Yli=i \ x i ~ Vi\ = 1- Let G = (Z d , Ed) be the infinite graph whose vertices are 
Z d , and whose edges, Ed, are the set of pairs (x, y) where x, y are neighbors. 

Note that the degree of each vertex is 2d. Set [m] = {1, . . . , m} for a positive integer m. Then [m] = 
[mi] x [7712] x ■ • • x [mj is a box containing # [m] = m\m2 ■ ■ ■ rrid points in Z d , where m = (mi, . . . , md) £ N d 
has positive integer coordinates. We view the finite lattice [m] as a graph G (m) = ([m] , _E([m])) where 
x, y £ [m] are neighbors, i.e. (x, y) £ -E (m) if and only if ||x — y|| = 1. Let T (m) be a d-dimensional torus 
on the finite lattice [m], where T (m) is obtained from the infinite graph G — (Z d ,Ed) by considering the 
quotient Z d / ((miZ) x • • • x (m^Z)). So T (m) is a 2d regular graph, and G (m) is a subgraph of T (m). 

A monomer-dimer tiling of [m] , also called a configuration <f> in [m] , consists of £ nonoverlapping dimers 
placed on the neighboring lattice points in [m]. Other lattice points in [m] are viewed as covered by 
monomers. Denote by C ([m] ,£) the space of monomer-dimer tilings <fi of [m] with £ dimers. Denote by 
Cp e r ( [m] , ^) all monomer-dimer tilings of T (m) with £ dimers. Clearly, C([m],£) C C per ([m],£). Let 
#C ([m] ,£) , #C per ([m] ,£) be the cardinalities of C ( [m] , £) , C por ( [m] , £) respectively. It was shown by 
Hammersley [19] 

W n log C ([mi] ,£i) \ogC per ([mi] ,4) 

X d (p) = hm — : = hm ^— = (2.1) 

m;^oo zfjz [m s ;J rrii^oo zffz [mi] 

where mi = (m^i, . . . ,mi i( j) £ N d , lim^oo m^k — 00, k £ [d], lim -m^] — P ^ [0>1]- This equality does 
not depend on a particular choice of sequences mi,£i,i £ N which satisfy the stated conditions. For more 
detailed analysis of the above limit and upper and lower bounds of Xd (p) see [17) . 

We now discuss the notion of p monomer-dimer entropy As (p) of a sequence S = {Si,i £ N} of finite sets 
in Z d . Assume that we have a sequence of increasing boxes [mi] , i £ N, where mi = (mi,(, . . . rrid.i) , i G N 
and for each k £ [d] the sequence m,-^, i £ N increases to infinity. Without loss of generality we may assume 
that Si C [mi]. Then Si = [mi] \E/j, for a sequence of subsets J7j C [mj , i £ N. Let = {([mi] ,Ui) ,i £ N}. 
It would be convenient to replace S by With Si we associate a subgraph Gi of T (mi) induced by the set 
of vertices Si — [mi] \C/j. (So Ui, Vj £ [mi] are neighbors if and only if they are neighbors in T(mj).) 
Let Q — {Gi,i £ N}. Then A5 (p) = Xy (p) is defined as equal to Xg (p) as in [TBI §2]- We now state this 
definition. 

Let [m] C N d , m = (mi, . . . , m d ) £ N d be a box in Z d . Assume that [/ C [m]. Denote by # ([m] \C7) the 
cardinality of the set [m] \U. Denote by C ([m] \U, £) all tilings of [m] \U with exactly £ dimers, 
when we view [m] as the graph G (m). (It is possible that #C ([m] \U, £) — 0). Then 

Aw (p) = 

r log#C([m t ]\[/ t ,^) 2^ 

hmsup — — — — provided hm — - — — — =p. (2.2) 

i^oo #([mi]\C/ l ) »^-cx) # ([m^ \Ui) 

Note that it is possible that Xu (p) = —00 for some p £ [0, 1]. We now state a conjecture, which seems to be 
very reasonable. 

Conjecture 2.1. Let d £ N, [mi] , mi = (m^i, . . . , nid,i) , i £ N be a sequence of increasing boxes [mi] , i £ N, 
where lim^oo TOfc,j = 00, fc = 1, . . . , d. Assume that Ui C [mj] ,i £ N anrf denote U = {([mj] , J7j) , i £ N}. 
T/ien 

Aw (p) < A d 0) for each p £ [0, 1] . (2.3) 
Proposition 2.1. Coniecture \2.1\ holds for d = 1. 



Proof. Clearly [m] \{7 consists of at most #=Z7 + 1 disjoint paths, Pi, . . . ,Pk- Let m! = # ([m] \{7) and 
view C? ([m']) obtained from Pi, . . . , Pj. by connecting by an edge the last vertex of Pi and the first vertex of 
Pi+i for i — 1, . . . , k — 1. Hence, each monomer-dimer tiling of [m] \U is a monomer-dimer tiling of [m'\. □ 

We assume that d > 2 and follow the discussion in 17 . Let m = (mi, . . . ,md) and denote m' = 
(mi, . . . , nid-i) ,N' = # [m'] , md = n. Then the (rf — l)-dimensional box of the form ([m'] ,j) is called the 
j th layer, for j = 1, . . . , n. So n — 1 corresponds to the (d — l)-dimensional case. Assume that n > 2 unless 
stated otherwise. We view the adjacency graph G (m) as composed of n layers of G (m') plus the edges 
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between the adjacent layers ([m'] ,j) and ([m'] , j + 1) for j = 1, . . . , n — 1. Let us consider one tiling of [m] 
with I dimers: <j) G C ( [m] , t) . Each layer ([m'] , j) has /ij horizontal, (transversal), dimers for j = 1, . . . , n. 
Let Vj be the number of vertical dimers connecting the levels ([m'] , j) and ([m'] , j + 1) for j = 1, . . . , n — 1. 
Let us assume that vq — v n — 0. We have the following conditions for n > 2. 

y)(/ij+«i)=4 Uj-i + Vj + 2hj < N' for i = l,...,n. (2.4) 
i=i 

More precisely, the location of Vj vertical dimers connecting the level j and j + 1 are ( V, ■ , j ') for a corresponding 
subset Vj C [m'] for j = 1, . . . ,n - 1. So Uj = #Vj. We assume that V = V n = 0. Then V^_i n Vj = for 
j = 1, . . . , n. Denote h = (hi, . . . , h n ), v = (v%, . . . , u„_i), «o = w„ = 0. Let C ([m] , h, v) C C ([m] , £) be all 
conhgurations of tilings of [m] , m = (mi, . . . , m^-i, n) with £ tiles, where h and v are the distributions of 
horizontal and vertical tiles. (It is possible that C ([m] , h, v) = 0.) Denote by C (m, £) the set of all possible 
(h, v) satisfying (|2.4j) and the corresponding above stated conditions. So 

C ([m] ,£) = U( h ,v)eC(m,*)C ([m] , h, v) , 

#C([m],*) = £ #C([m],h,v). (2.5) 

(h,v)eC(m,£) 

We claim that 



Clearly /^ G 



^ 

' 9 



#C(m,l)<(N'f n . (2.6) 
, G [0, N'j. Hence the number of choices of hi, is at most (N' + 1) which is not 



greater than (N') 2 . So (|2.6|) is a big overestimate, since we ignored ()2.4|) . 

We now show that to compute \d(p) we need only to consider one sequence hj,Vi,fc G N of dimer 
distributions. 

Lemma 2.1. Let d > 2 be an integer. Assume that m^ G Z d , ^, i G N satisfies the conditions for (|2.ip . XTien 
£/iere exists a sequence 

h, = (/l M , . . . , ft mMii ) G Z™*' d , v, - . . . , U^-ij.i) G Zt" 1 ^- 1 ), 

«»,o = Wj.nn.i = 0, ^ hj,i + Vj,j =£j, i € N. (2.7) 
i=i 

suc/i £/ia£ 

lim io g# c(K] ,h,,v <) = 
i^oc # [niij 

Proof. Let fhi,ViJ , i G N be an allowable sequence satisfying (|2.7p . The equalities (|2.1I) and (|2.5[) yield 
the inequality 

log#C ([mi] ,hj,Vi) 
limsup -^-^ '- < X d (p) . (2.9) 

Choose 

(h», V{) = argument of max j#C (liOi] ,&i,v^ , (hi,ViJ G C (mi,4)j ,i G N (2-10) 
i.e. a value of fhj,Vj^ for which maximum is achieved. Combine (|2.5[) and (|2.6p to deduce that 

#C ([mi] , Q < (# [m',]) 2mi ' d #C([mi] , h 4 , Vi ) , i G N. (2.11) 



Observe next that 



hm M#KIp! = um 2iog#K0 = o (2 . 12) 
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(l2"TTD and (l2"T2l yield 



Combine this inequality with (|2.9[) to deduce the lemma. 



log #C ([mi] 

liminf — = '- > X d (p) . (2.13) 



3. A Lower Estimate of \d(p) 
The main result of this section is: 
Theorem 3.1. Let d > 2 be an integer and p G [0, 1]. Then 

A d (p) > (Vi (g) + (1 - g) log (1 - «,) - log + (3.1) 

i - ^y^J lo s - ^~y^) ~ ^ ~~ ^ log ^ ~~ for each g G ' 

In particular (|1.12[) holds. 

We prove the theorem using a number of known results and straightforward lemmas. Recall that A^ (p) > 
Ad-i (p). Hence to show (13. 1|) we need to consider the case q < p. Next we recall the inequality 

V#[ m ]/ #H 

This follows by considering boxes [jm] tiled with j d l dimers for j = 1, . . . , 117j . 

Fix a rational number q G (0, 1) and e > 0. Then there exists m' G N d_1 ,^' G N such that 

21' w , v log C ([my) 

9:= #RI' Arf - l(g):= #m >^fa)- e - ( 3 - 3 ) 

See [15]. Let 

or 

p = jT-p — where 2/ < 2L < # [m'] . (3.4) 
We will show for these choices of p and q the following inequality holds. 

A d (p) > (a^ (g) + (1 - g) log (1 - g) - log + (3.5) 



^)log(l-^)-(l- P )bg(l-p) 

Since A& (x) is continuous for each x G [0, 1] and each k €N [T5], this will imply (|3.1[) . 
Lemma 3.1. Let IA be the set of all subsets of [m'] of cardinality L — £' . Then 

#C(K] , O = /o X! # c (M W> . # w = ( jZ e ) ( 3 - 6 ) 

N' - 21' 



where N' = # [m'] , / = 



Proof. Take a tiling of [m'] by I' dimers. The monomers of this tiling are located on the subset V of 
[m']. Clearly #V = N' — 21'. This tiling appears exactly in all C ([m'] \U,£'), where U C V. The number 

N' -21' 



of choices of these U is ^ ^ ^ J . The formula for jftA is clear. We have implicitly used the fact that 

L-l' <N' - 21'. " □ 

We now consider the d-dimensional box [m], where m = (m',n) and n > 2. So [m] is viewed to have n 
horizontal levels consisting of [m']. Let h„ = (£',...,£') G N",v„ = (L — £',..., L — £') G N" _1 . For each 
(7eM we denote by C ([m] , h„, v„, U) a subset of tilings in C ([m] ,h„,v n ) such that no dimer (vertical or 
horizontal) covers a spot in U on the first level ([m'] ,1). 
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Lemma 3.2. Let n > 2, m = (m',n), h„ = (£',...,£') £ W\ v„ = (L —£',..., L — £') £ N n . Then 

#C([m],h„,v n ) = /^#C([m],h„,v„,C/), where f=( N ' ~ ^ + l> ) \ (3.7) 



Proof. Take a tiling <f> g C ([m] , h„, v„). The monomers of this tiling, located on the first level ([m'] ,1), 
form a subset W of [m']. Clearly #1U = N 1 — (L + £'). This tiling appears exactly in all C ([m] , h„, v n , U), 

( N' — (L + £') \ 

where U C W. The number of choices of these U is I ^ ^ J . □ 
Lemma 3.3. For n G N 

#C ([(m', 2n)] , h 2n , v 2n ) = £ #C ([(m', n)] , h„, v„, C/) 2 , (3.8) 

ueu 

where C([(m',l)] ,hi,vi,C/) := C ([m'] \f7, £')• 

Proof. Observe that any tiling in <j> £ C ([( m 'i 2n)] , h 2 „, v 2n ) can be obtained from exactly two tilings of 
C ([(m', n)\ , h n , v„, U), where U is the location of L — £' vertical tiles from level n to level n + 1 in 4>. □ 
Using Lemmas 13. II - 13.31 (13.31) and the Cauchy-Schwarz inequality we obtain: 



Corollary 3.1. 

e 2 



e 2AT'A' d _ 1 (g) / N , 

#C([(m',2)],h 2 ,v 2 )> — — 5 , M = m=[^_ l , , (3.9) 



ji/^/r/' / r> \i 1 x . (#C([(m',n)] ,h„,v„)) 2 

#C([(m ,2n)] ,h 2n ,v 2 „) > , for n > 2. (3.10) 

Let s€N and consider C ([(m', 2 s )] , h 2 s , v 2 « ). Use the above corollary to deduce 

#C([(m',2-)],h 2S; v 2s )>Af/ 2 l Mfh j • (3.11) 

2(2 S L— ( L— ^ 

Note that each configuration in C ([(m', 2 s )] , h 2 s , v 2 s ) has 2 S L — [L — £') dimers. Let p s = — — — — . So 
limg^oo p s = p. (|3.2p yields that 

AiWjM (m^^ (312) 

Combine the above inequality with (|3.11l) and let s — > 00 to deduce 

A d (p)>A' d _ l(9 )-^^- (3.13) 
Let N' —> 00, i.e. m' — > 00, while keeping the values of g and p fixed, to deduce (13. 5p . 

4. An Upper Estimate of Xd{p) 

In this section we use the notation from Section [2j We will study a sequence of cubes and use Lemma 12.11 
so that with each of these cubes, say m^, we need consider only a single and Vj. We turn to considering 
a single such cube, m, suppressing the index i since we spend much time with this single cube. 

Recalling notation from Section [3J 

m = (mi, . . . , md) , here all m, = n (4-1) 
m' = (m 1 ,...,TO d _i) (4.2) 

h=(h u ...,h n ) (4.3) 
v =(«i,..., «„_i) (4.4) 
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hi is the number of horizontal dimers in the i th layer, and the number of vertical dimers connecting the 
i th and (i + l) th layers. Vi is the location in m' of the vertical dimers connecting the i th and (i + l) th layer. 
We set as before 

N' = #m' = n d - x (4.5) 
and take I to be the total number of dimers. We then have the following. 

V,nV l+1 =0 i = l,...,n-2 (4.6) 

v t +v t+l <N' i = l,..., n -2 (4.7) 

2hi + + Vi < N' % = %..., n-2 (4.8) 

(There are slight obvious modifications of (|4.6p - (|4.8I) for i = l,n — 1, n.) 

For our given v we count the possible choices of V\, . . . , V n -\ such that = Vi for i = 1, . . . , n — 1. We 

can choose Vi in ( | ways, Vo in | Vl | , etc. Thus we arrive at 

)n ("*--•) ( «o) 

for the total number of choices. Using this enumeration, and recalling Lemma 12.11 we will get an upper 
bound for Ad (p) as the limit as n —> oo of the maximum possible value of 



flog a (n d -\ v) + log J] C ([m'] \V-_x U V^h^ 



(4.11) 



We define 



There follows 



2/ 

P=~ d (4-12) 
n 

^ = ^TT' *< = ( 4 - 13 ) 

^ n 1 n 

u=-y Ut , i=-Vf, (4.14) 
n * — ' n z — ' 

i=l i=l 

~ _ 2hi 

U n^(l-\(u l+ u^)) (4 - 15) 

u + t = p (4.16) 

u is the density of vertical dimers and t the density of horizontal dimers. tj is the fraction of the vertices in 
[m'] \ (Vi-i U Vi) covered by dimers. 

Using Conjecture 12. II we may effectively assert 

logCtKlU^uv^A) < 

Hence 

ipgnsu c(Ki \ w- u y.) a) . ^ i - gag Xji (t) (41S) 

i=l 

Recall that Ad-i (p) is a concave function. (See [16] for a sharper statement.) Since ~ Y^7=i (l 
1 — it it follows that 



^ (tj) < (1 — u) A d _i (t 1 ) , t' = Y^— 2 — U. (4.19) 

1 v 7 



i=l 



s 
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In view of (4.15) we obtain 



1 - 2hi 



t' = Y-, 2 v = -, L—Y2h l = -— -. (4.20) 

^ (l-u)n n rf-i ( x (l-w)n d ^ (1-u) v ; 



We combine with the previous inequalities to deduce that 



(l-u) 

To bound the first term in (|4.11l) it is straightforward to show the maximum is achieved with all Vi equal, 
v t = This yields 

limsup ^^lA <(!_«) log (l - |) - | log | - (1 - „) log (1 - „) . (4.22) 

In conclusion we derive the following. 

A d (p) < max { (l - |) log (l - |) - | log | - (1 - u) log (1 - «) + (1 - u ) A<j-i (^-^j } (4.23) 

Because the expected value of u is § , by the symmetry of the cube, one can show one can restrict the range 
of it in taking the maximum in (|4.23|) to either [0, ?] or This leads to two formulas for an upper 

bound. Unless these are equal, one gets a better bound restricting to one or the other of these ranges. 

One can employ the argument behind Lemma 13. II in the construction of this section, obtaining an alternate 
upper bound: 



\d (p) < max |-- log - - (l-p) log(l-p) + [1 -p+ -J log (l 



u 

" + 2 



Again the range of u may be restricted to either [0, ?] or [§,l]- We outline this alternate development, 
closely following the previous route, and using the same notation. 

We let W r (V r ) be the number of configurations in the portion of the lattice living in layers r + l,r + 
2, . . . , n (summing over both vertical and horizontal dimers). This will involve summing over V r +\, . . . , V n —i 
compatible with # (Vi) = Vi,Vi D Vi+\ = 0, and hi horizontal dimers in layer i for i = r + l,...,n. 
V r , h r+ \, . . . , h n , Ur+ii ■ • • j v n-i are fixed. Then we have 

W n -x (K_i) < C ([m'\ \V n -x,h n ) 

< ^(l-Un-^,-^) = (4 25) 

with 

£ = 2/tt (4.26) 

We find the bounds bi inductively. Assume we have bounds 6 r +i, 6 r +2, • • • , &n— l- We seek 6 r . 

W r (K)< ^ C([m']\7 r UU h r+1 ) 6 r+1 (4.27) 

Vr+i 
#(K+l) = Wr+i 

r\V r = 

By the argument of Lemma 13.11 

< C ([ml \F r , h r+1 ) fe r+1 f ^ " 2 ^ +1 ~ Wr ) (4.28) 
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As before 



x d (p) < i • n 

>-E;(i-?)^(M 



which by a counting argument 



and 



< 



1 — - 

2 



(4.30) 
(4.31) 



(4.32) 



A' - 2h r+1 - v r 

V r+ l 

The maximum in II occurs where all Vi are equal and all h{ equal (at least for i = 2, . . . ,n — 2, just as good). 
This leads to (H1M1) . 

We now study the relation between ([TT3"j) and ([TT4]) (or (l4~23"]l and (|4^4]) 1. We take \ d (p) to be given 
by its approximate form 



Ad (p) = ^P log (2d) - | logp - (1 - p) log (1 - p) - | 



(4.33) 



and substitute this into 



(1 - 1) * (1 - 1) - 1 * (|) - (1 - -) lo S (l - .) + (1 - .) A,-, (^) } - 

(l -P + |) log (l + |) - | log (|) - (1 - P) logfl -p) + (l - I) Ad_i (f^|) } (4.34) 



Surprisingly, by a simple computation this comes out to be 



i (p - u) log 



1 _ « 

1-tt 



(4.35) 



In the region u < p this is positive. This shows that the expression in braces in (|1.13[) is pointwise (in u) 
greater than the expression in braces in (|1.14j) . Thus (|1.14j) is a better upper bound than (|1.13[) for (p) as 
given by (|Q5j> . 

For < u < p/d and d > 2 we have 



- (p - u) log 



I - a 

2_ 

l-U 



> c\pu 



We let the actual Ad (p) be related to the approximate form, (14.33[) . as 

A d (p) = Ad (p) + E d (p) 



We assume for d > do one has 



and 



\E d (p)\ <C2~ d 



dp 



E d (p) 



<c 3 - 



(4.36) 

(4.37) 
(4.38) 

(4.39) 



Then there is a d such that if d > d then (I1.14j) provides a better upper bound than (|1.13l) . That is if one 
substitutes (for d > d) (|4.3T[) into the right side of (11.131) one gets a larger value than if one substitutes into 
the right side of (|1.14[) . One has used the bound 



(1 - u) £d-i 



p — u 
l-u 



E d -\ 



p — u 
1 — — 

2 



< 



CiU |^d-i (p*)| + c 5 pu \E' d _ 1 (p**)\ 



(4.40) 
(4.41) 



Provided (|4.38|) and (|4.39|) are satisfied by E d (p) it is easy to deduce, using (|4.36|) and (|4.40|) . that (|1.14|) is 
a better bound than (|1.13[) for d large enough. One looks at the u* that maximizes the expression in braces 
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in (|1.14[) . By (|4.36|) and (|4.40l) . for large enough d, the expression in braces in (|1.13l) evaluated at u* is larger 
than the expression in braces in (| 1 . 14[) at u*. Thus the bound (|1.14[) beats the bound (|1.13[) . 

5. The Asymptotic Expansion 

This section largely parallels the development in with modifications appropriate to the monomer- 
dimer problem. We do not assume familiarity with |llj . 

We work with a periodic cubical lattice, A, of even edge length. We denote the volume, or equivalently 
the number of vertices, by either V oy N . A fraction p of the vertices are covered by dimers, thus pN/2 
dimers are used in each covering. We will be using the terms 'p-tiling' and 'p-cover' to describe coverings of 
a fraction p of the vertices. 

In d dimensions there are d 'kinds' of dimers, each kind oriented in one of the d lattice directions. Each 
kind of dimer may be 'located' at some place on the lattice. We generalize this situation as follows. A 'located 
tile' is a two element subset of the lattice. A 'tile' is an equivalence class of located tiles, with equivalence 
given by setting two subsets equivalent if one is a translation of the other. (A 'tile' is a generalization of a 
'kind of dimer'.) Notice that tiles need not be connected. 

We consider tiles with a 'weighting', a function on tiles. We normalize the weightings we consider, by 
requiring, if g is the weighting function, 

= 1/2, (5.1) 

t 

the sum over all tiles, t. We let / be the weighting function given by 

[ to if t is a dimer, 
f(t) = < (5.2) 
otherwise, 



that clearly satisfies the normalization condition (|5.1I) . 
A 'p-tiling' Ti of A is a set of two element subsets of A, 



Ti = 



| s 1j s 2i ■ ■ ■ ) S pjV/2} ) (5-3) 



where the s\ are disjoint, they are located tiles. 

We now realize the sum over all possible dimer p-covers of A, the goal of our study, as 

{2d) pN/2 Z (5.4) 

with 

z =/Z n ( 5 - 5 ) 

where the sum is over all p-tilings of A; the product over the /'s selects those tilings in which all the tiles 
employed are dimers. The bar over a subset indicates the equivalence class of the subset as defined previously. 

We let fo be a constant function on tiles with value , to satisfy the normalization condition f|5. 1[) . We 
write 

/ = fo + (f- fo) =fo + v, (5.6) 
/, /o, and v all functions on tiles. Z becomes 

Z = E II (/o + «(««)). (5-7) 

Z = Z + Z x + Z 2 + • • • (5.8) 

having expanded Z in powers of v. 
We note 

e NX *<*> = (2d) Np/2 Z. (5.9) 
Here (p) is understood to be a function of N, the usual Ad (p) the infinite volume limit 

A d (p) = |log(2d) + -ilogZ. (5.10) 
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If one replaces Z by Zq, the mean field approximation in a natural nomenclature, one gets taking the infinite 
volume limit 

\ d (p) £ | log (2d) - | logp - (1 - p) log (1 - p) - | 

= flog (2d) + lim ±-\ogZ (5.11) 

by an easy calculation. 

We return to (15.81) and introduce some convenient notations: 



Z = Z a Z*, (5.12) 

Z* = 1 + Z{ + Z* 2 + ■ ■ • , (5.13) 

Z* = Zi/Z . (5.14) 

There is a natural factorization of Zi into a contribution from the factors of v in (|5.7p which we call Z* and 
the factors of fa in (|5.7p which we call f3 (N, i) Zq so that 

Z*=fi{N,i)Z* (5.15) 

with 

P(N,jN)~e NH(p < j) (5.16) 

where 

H (p, j) = (1 - 2j) log (1 - 2j) + | logp - (| - i) log (p - 2j) . (5.17) 

Equations (|5.16|) and (|5.17|) follow from a short computation, always working in the large N limit. 
We let Z* be Z* with (N,i) replaced by 1, 

Z* = 1 + Z\ + Z; + ■ ■ ■ , (5.18) 

where a detailed specification of Z* is given by 

2 i=\ e n ( 5 - 19 ) 

Si,S 2 ,.-.,Si a=1 

disjoint 

Now referring to [3] we may write a cluster expansion for Z* 

logZ* = ^ij s , (5.20) 

J s = ^2 v (si)---v(s s )tfj' c (s 1 ,s 2 ,.-.,s s ), (5.21) 

si,s 2 ,...,s a 

where we may identify Eq. (2.5a) of [3] with Z* , and ([5301) . (f5T2l]) with Eq. (2.7) of j3J. The located 
tiles in the sum of (|5 . 2 1 1) are forced to overlap so they cannot be divided into two disjoint sets. ip' c is a 
numerical factor depending on the overlap pattern. To make our computations mathematically rigorous it 
will be necessary to study the convergence properties of sums such as in (|5.20l) . At present this appears very 
difficult, and we by no means see yet a clear route to a proof. It is a challenging problem for the mathematical 
physicist. 

It is easy to show J\ = and it is proven in [9] that 

J s = —, rn — h---H ; — t- (5.22) 

d r d r+1 d 3 - 1 

with r > s/2. We also find it convenient to define 

NJ t = (1/i!) Ji. (5.23) 

From (|5TT5|) , (j5~T5|) . ([5TS]) . (j5~2^| one gets 

Z*= E K^-)^-^ 1 ^- (5-24) 

Qi,...,a s + i 
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We approximate the sum in (|5.24[) by its largest term, in the limit N — > oo. If all the J's are positive this 
is a reasonable way to extract the dominant asymptotic limit. In [10] an argument is given that our results 
will hold even if some of the J's are negative. 

Most important to observe is that the J's of this paper are the same as the J's of the dimer problem! The 
Ji were computed through J 6 as follows, see [TT] : 



Ji = 0, 




(5.25) 


J2 = \\' 




(5.26) 


Js--- 

3 ~~ 12 d 2 ' 




(5.27) 


- 3 1 
M ~ 32 d 2 " 


3 1 
+ 64 W 


(5.28) 


7 - 1 1 


3 1 

" 80 d*' 


(5.29) 


7 - 7 1 
6 ~ 48 d 3 


5 1 11 
64C7 1 ~ g^5- 


(5.30) 



The computations were done in integral arithmetic using Maple. Because of (|5.22[) . to compute the Ji 
through Jq for all d, it is sufficient to compute these Ji for one, two, and three dimensions. The most 
complicated computation of these involved placing down six dimers in three dimensions, in all possible 
ways where overlaps make it impossible to disconnect into two disjoint subsets of dimers. One is using an 
exponential time algorithm, Js's computation took 3 seconds and Je's two weeks. To compute the 1/d 4 terms 
in (|1.2[) would require knowledge of J7 and J% also, we think unless a new method is found one will never do 
this computation! 

We return to computing the largest term in (I5.24[) . We differentiate with respect to the a's in (|5.24[) . using 
(|5.16[) and (|5.17[) to deal with the factor of (3. Finally scaling Oj — y jron we find 

loga fc = log J fc + (Pi ^2 ia *) = log ^ k + Fk (5.31) 

and 

Z* ~ e N {-'E a i F i+'£ J i eFi + H (p^ ia i)}, (5.32) 

From these two equations, and (|5.10p . (|5.1ip . (|5.12[) , one can find (p) as a formal power series in the J's, 
and then a formal power series in 1/d using the expressions (15.25)) through (|5.30[) for the J's. This algebra 
we did using a 40- line Maple program, resulting in (|1.5p - (|1.8p . 

In [12] the first named author shows that if the terms in the formal expansion (jl.5p are rearranged as a 
power series in p, then for sufficiently small p this series converges. 

6. The LAMC and Other Conclusions 

For any infinite graph G = (V, E) on a countable number of vertices V, where each vertex V has at most 
degree r, one can define the notion of the p-entropy \g(p), similarly to the definition of the entropy Xu(p) 
defined in §2. The Lower Asymptotic Matching Conjecture for infinite regular r-bipartite graphs claims that 
\g(p) is bounded below by a universal function u) r (p), see equations (1.3) and (1.4) of [16 j. uj r {p) can be 
considered as the p-entropy of a random r- regular infinite bipartite graph. The formula for uJ2d(p) is given 
by the right-hand side of (|6.ip . We first prove the LAMC for a rectangular lattice, 



Theorem 6.1. The LAMC is true for a rectangular lattice. That is 
1 



2 



A d (p)>- plog(2d)-plogp-2(l-p)log(l-p) + (2d-p)log l-£ (6.1) 



P 



2d 



Proof. We prove this by induction. It is true for d = 1 by [TB] eq (1.1) or [T7] § 4. In the inductive step 
we are given 

1 " 



Ad-i (p) > 2 



plog (2 (d - 1)) - plogp - 2 (1 -p) log (1 - p) + (2 (d - 1) - p) log ( 1 - '" 



2(d-l), 



(6.2) 
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We substitute this into (the right side of) (|1.12[) . proven in Theorem 13. 11 and arrive by a simple computation 
at exactly (|6.1[) . 

We turn to a study of the relation between our upper and lower recursive bounds and the asymptotic 
expansion. We assume that A^ (p) is given by an asymptotic expansion in inverse powers of d, whose first 
few terms are as in (| 1 . 10|) . 

A d (p) ~ \v log (2d) - | log (p) - (1 - p) log (1 - p) - | + ^ + hp ^f +... (6.3) 

Theorem 6.2. 

A) a>\ 

B) Ifa = i, then (-2 + 2bp + c) > 

C) Ifa=\ andc = 2, then b>0 

D) Then values ofa,b,c as given in p. 61) and (jl.7|) ensure that the asymptotic expansion in (| 1 . 5|) satisfies 
the implications of the lower recursive inequality (jl.121) . 

Proof. If we substitute (|6.3[) into the recursive inequality (|1.12p . then the value of the right side minus 
the left side is given in (II. lip which we rewrite 

-i<- 1+8 °>$-^i- 2+2f * + 4 +0 (i) <"> 

The requirement that (|6.4I) is asymptotically < yields parts A, B, and C of the theorem. The values of 
a, b, c as given in (|1.6[) and (|1.7[) ensure (|6.4|) is asymptotically negative. 

In short our asymptotic expansion for Ad (p), (|1.5|) . satisfies the lower recursive inequality. 

The upper recursive inequalities we have (|1.13[) . (|1.14p are not strong enough to put conditions on the 
coefficients of the asymptotic series, such as parts A,B, and C of Theorem 16.21 

7. Numerical tests for d = 1, 2, 3 

In the last section we saw our expansion satisfied the implications of the lower recursive bound asymptot- 
ically. Here we study this expansion and our bounds for small d, d = 1, 2, 3. 

We find it convenient to rearrange our expansion for A<j(p) as a power series in p, perhaps this is always a 
better form to work with. 

X d (p) ~ 1 (pln(2d) - plnp - 2(1 - p) ln(l - p) -p) + Md)p k (7.1) 

k=2 

where we see from (1 1 . 5|) - Q1.8P that 

a 2 (d) = H (7.2) 

as(d) = 11 (7.3) 

a4(d) = 32^-192^ (7 - 4) 

using the fact mentioned after (|1.5p that c/c (p) is built up of powers p s with k < s < 2k. Knowing Ji for 
i < 6 we can also compute: 

1 1 39 1 

a5(d) = l6^ 640^ (? - 5) 

a6(d)= 24^-32rfJ- 1920^ (7 ' 6) 
7.1 d=l 

Recall the exact formula for Ai(p), e.g. (TTl end of §4] 

Ai(p) = (l - |) log (l - |) - | log (|) - (1 -p)log(l -p). (7.7) 
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Expand log(l-f) in. power series in p to deduce that 

1 00 k 

A x (p) = -plog (2) - | log (p) - (1 - p) log (l-p)-f +E {k -l)k2« - ( 7 - 8 ) 

It is easy to check the values of afc(l) are correctly given to match (|7.8[) for fc = 2, . . . ,6. (We have not 
tried to prove this to all orders, the Ji for d = 1 can all be computed, but it probably can be done.) This is 
a remarkable validation of our procedure to compute A<j(p). 



7.2 The Expansion for A 2 (l) and A 3 (l) 

In eq. (I7.1[) we consider putting the upper limit in the sum to be n. Then for n — 2, . . . , 6 we get a 
sequence of "approximations" to Xd(p)- Setting d — 2, p — 1 we get the sequence 

.2556, .2609, .2654, .2694, .2724 (7.9) 

and setting d = 3, p = 1 we get 

.4375, .4399, .4424, .4439, .4450. (7.10) 

The exact value for A 2 (l) is given as 

A 2 (l) = .2915... (7.11) 

by [14] and [23], and A 3 (l) satisfies 

.440075 < A 3 (l) < .457547 (7.12) 

by @] and [TBJ. We may well believe the sequence (|7.9I) extended converges to A 2 (l) and (|7.10|) extended 
converges to A 3 (l). We now believe (|7.1[) converges for all physical p, not just small p as proved in [T2"] ! 



7.3 The Expansion and Bounds for A 2 (p) 

We present a table of values for A 2 (p), its bounds and approximations. The first column gives the value 
of p. The second column, the expansion value, is the result of keeping in (|7.ip the dfc(2) from k = 2 to 
k = 6, the known values. The third column is the rigorous lower bound (|1.12[) . The fourth column presents 
the "exact" value for A 2 (p), actually a very good approximation from [17] and pQ. Looking at the table 
on page 654 of Baxter's paper pQ, a row with the three column entries cl,c2,c3 corresponds to p = 2c 3 , 
A 2 (p) = In (c 2 ) — (1 — p) In (ci). The values of p thus in this table led to the choices for p in our table. The 
final column is the upper bound from (|1.14j) . actually rigorous for d = 2. 



p 


exp 


lb 


exact 


ub 

















.14870 


.30887 


.30887 


.30887 


.31030 


.26030 


.45283 


.45281 


.45284 


.45734 


.50426 


.63492 


.63449 


.63495 


.65274 


.77053 


.62983 


.62678 


.63086 


.67319 


1 


.27236 


.26162 


.29156 


.34657 



We are struck by how much better the lower bound is than the upper bound, and how good an approximation 
is the expansion, keeping the terms we know, for p < .5. 



8. Future Directions 

We leave to future research three compelling problems: prove that the "asymptotic" expansion expressed 
as a power series in p converges for all physical p and d, find better upper bounds pure or recursive, prove 
Conjecture 12.11 
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